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Abstract 

Motivated by holography we explore higher derivative corrections 
to four-dimensional Anti-de Sitter (AdS) gravity. We point out that 
in such a theory the variational problem is generically not well-posed 
given only a boundary condition for the metric. However, when one 
evaluates the higher derivative terms perturbatively on a leading order 
Einstein solution, the equations of motion are always second order and 
therefore the variational problem indeed requires only a boundary con- 
dition for the metric. The equations of motion required to compute 
the spectrum around the corrected background are still generically 
higher order, with the additional boundary conditions being associ- 
ated with new operators in the dual conformal field theory. We discuss 
which higher derivative curvature invariants are expected to arise in 
the four-dimensional action from a top-down perspective and compute 
the corrections to planar AdS black holes and to the spectrum around 
AdS in various cases. Requiring that the dual theory is unitary strongly 
constrains the higher derivative terms in the action, as the operators 
associated with the extra boundary conditions generically have com- 
plex conformal dimensions and non-positive norms. 
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1 Introduction 



In this paper we will explore higher derivative corrections to gravity theories 
in (3 + l)-dimensions with negative cosmological constant. Our main moti- 
vation for looking at higher derivative corrections to four dimensional AdS 
black holes is in the context of holography and, in particular, applied holog- 
raphy, AdS/CMT, where many of the systems of interest are modelled by 
four dimensional bulk spacetimes. The addition of higher derivative terms 
allows us to probe the dual physics as one moves away from infinite N and 
infinite 't Hooft coupling. 

Finite N effects can change the physics qualitatively. For example, let us 
consider holographic superconductors, a subject which has been extensively 
studied in recent years, initiated by [1], [2] and [3]. Working with classi- 
cal gravity there is an apparent violation of the Coleman-Mermin- Wagner 
theorem. This well-known theorem states that, for system in two spatial 
dimensions, we cannot have continuous symmetry breaking in systems at 
finite temperature and hence the formation of a symmetric breaking con- 
densate is forbidden. However, holographic superfluids have been found in 
(3 + 1) bulk dimensions, in which a symmetry breaking operator in the dual 
(2 + 1) dimensional CFT acquires an expectation value. As explored in [4], 
this is an infinite N effect and at finite N quantum effects in the bulk indeed 
ensure that the symmetry breaking operator does not have a well defined 
expectation value, in accordance with the expected field theory behaviour 
[5]. 

One does not see a qualitative finite N effect such as the restoration 
of the Coleman-Mermin Wagner theorem by evaluating higher derivative 
corrections on the leading order gravity solution but rather by exploring 
quantum effects in the bulk. Evaluating higher derivative corrections rather 
shifts the saddle point and allows one to compute corrections to thermo- 
dynamic quantities, transport coefficients and so on. In the context of five 
bulk dimensions, a considerable effort has been put into investigating higher 
derivative corrections and exploring the effects on the ratio of the shear 
viscosity rj to the entropy density s, see for example [6, 7, 8, 9]. 

In particular, [9] used Gauss-Bonnet curvature corrections and initiated a 
bottom up exploration of the constraints on the higher derivative corrections 
imposed by unitarity of the dual CFT. Working with the Gauss Bonnet term 
is particularly convenient because the corrections to AdS planar black holes 
are known analytically for any value of the Gauss Bonnet coupling constant, 
see [10] and also [11, 12, 13]. Note that an effect of the Gauss-Bonnet term 
relevant to the superfluids mentioned above was discussed in [14], where it 
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was found that addition of the higher curvature terms makes condensation 
to a superfluid phase more difficult. 

The Gauss-Bonnet terms, and corresponding corrected AdS black holes, 
are a useful way to go beyond classical gravity in bulk dimensions higher 
than four. However, such terms are trivial in four bulk dimensions, in the 
sense that an Einstein metric is uncorrected and therefore one needs to 
include higher order curvature invariants to obtain non trivial corrections 1 . 
An alternative possibility is to couple Einstein gravity to a dilaton in four 
dimensions because Gauss-Bonnet type corrections to diatonic black holes 
are then non trivial, see for example [16], but this does not address the 
question of how AdS black holes with no dilaton are corrected. 

Apart from AdS/CMT motivations mentioned above, for which dilatonic 
AdS black holes may indeed already capture many relevant features, see for 
example [17], there are a number of other important motivations in exploring 
higher derivative corrections to Einstein gravity with a negative cosmological 
constant. The first is in understanding the AdS/CFT correspondence when 
the dual theory is on an S 3 . In recent years there has been considerable 
progress in understanding dual (supersymmetric) 3d CFTs, following the 
works of BLG [18] and ABJM [19], and localisation techniques have been 
used to compute free energies of the dual theories placed on an S 3 . Taking 
the limit of large ./V and large 't Hooft coupling, the free energies have been 
matched to the onshell renormalised action of AdS^ with an S 3 boundary 
in Einstein gravity [20]. Localisation techniques also allow us to access the 
subleading terms in the free energy which should be compared to the effects 
of higher derivative terms evaluated on the bulk AdS±. Comparing these 
subleading terms with the gravity results we develop here can be used to 
test the correspondence and indeed restrict which higher derivative terms 
can arise in the four dimensional bulk action. 

The second motivation in exploring higher derivative terms in four di- 
mensions is in the context of understanding the holographic dictionary. One 
of the main points of this paper is that the addition of higher derivative 
terms generically involves additional data being required for the variational 
problem to be well-defined. In the context of holography, the additional data 
corresponds to a new operator in the dual CFT, in addition to the stress 
energy tensor which is dual to the bulk metric. For generic higher derivative 
terms added to the action the dual operator has complex dimension and/or 
negative norm, reflecting the fact that the corrected added violates unitar- 

1 A review of higher order gravity theories and their black hole solutions may be found 
in [15]. 
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ity. This analysis provides a very direct probe of the unitarity properties 
arising from the higher derivative terms. 

Historically the main context in which higher curvature corrections to 
four-dimensional gravity has been studied is as a toy model for a quantum 
theory of gravity. In this context the key problem is that adding higher cur- 
vature corrections adds higher-order time derivatives to the theory and con- 
sequently ghosts. Recently there has been considerable interest in so-called 
critical gravity theories, in which ghostlike modes appear to be absent, in 
both three [21] and four [22, 23] bulk dimensions. The four dimensional story 
that we develop here is the exact analogue of the discussions in [24, 25] for 
topologically massive gravity in three dimensions [26, 27]: the higher deriva- 
tive terms in TMG were shown to be associated with a new operator in the 
dual two dimensional CFT. In TMG, regardless of the value of the coupling 
of the higher derivative term a violation of unitarity was found in the dual 
field theory, either by a complex operator dimension or by an operator whose 
two point function was non-positive. Note that this violation of unitary per- 
sisted even at the so-called critical point, where the new operator together 
with the stress energy tensor were non-diagonalizable. In this paper we 
will show that analogous problems are found in the four-dimensional higher 
derivative theories. 

Given that higher derivative terms generically give rise to new boundary 
conditions and hence dual CFT operators, whose properties are not consis- 
tent with unitary, one may ask how this observation can be consistent with 
the fact that top down models arising from string theory are necessarily uni- 
tary. To understand this point, one should first note that in the context of 
string theory and holography the higher curvature terms are always viewed 
as an infinite series of small corrections to the leading order effective action. 
The action with higher derivative terms is not quantized, which makes the 
issue of ghostlike modes moot. In other words, the effective action takes the 
form 

where A is the cosmological constant; R n denotes schematically an n-th 
order invariant 2 ; a n is a dimensionless numerical constant and l p denotes 
the effective Planck length. The effective Newton constant in the Einstein 
theory is k 2 = 8irG. 

One is by assumption working in a regime where l p is small and therefore 
the corrections should be treated perturbatively. Suppose <7( ) is a solution 

derivatives of the curvature can also arise but will not be considered here. 
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of the Einstein theory, namely 

0^(9(0)) = R^{9(0)) ~ ^R(9(o))9(o)uv + As>( )^ (2) 

Then the corresponding solution of the corrected theory can be expressed 
as a perturbative series 

9 = 9(0) + lp +l 9{n) + ■■■ (3) 

with 

G^(9(n)) = "On ^£7 (0(0)) (4) 

and so on. 

One should emphasize at this point the conceptual difference between 
evaluating the higher derivative terms on the lowest order solution and treat- 
ing the higher derivative term non perturbatively. In the former case, the 
equations for all the metric corrections gi n ) are second order inhomogeneous 
differential equations, rather than higher order differential equations. Since 
the equations are second order, the only boundary data that needs to be 
supplied for the variational problem to be well-defined is the metric. When 
one is considering the higher derivative terms evaluated on the lowest order 
solution, an analogue of the Gibbons-Hawking- York [28, 29] term in the ac- 
tion can always be defined such that the variational problem is well-defined 
for a Dirichlet condition on the metric. 

By contrast, as we will explore in sections 2 and 4, whenever the higher 
derivative terms are treated non-perturbatively or when we consider the 
spectrum around a given corrected background, the resulting equations of 
motion are generically higher order 3 . This means that additional boundary 
data needs to be supplied. In the context of holography one can understand 
the additional data as corresponding to additional dual operators in the 
field theory, beyond the stress energy tensor. The variational problem in 
such contexts will be well-defined only if one supplies additional information 
together with the Dirichlet condition on the metric; the actual information 
which is needed depends on which higher derivative terms are added. 

Thus, given a background which solves the supergravity equations at 
leading order, the variational problem will be well-defined when one com- 
putes the corrections to this solution without specifying additional data. 
However, when one looks at the spectrum around this background, the 

3 The Lovelock theories [30] are a well-known counterexample in which the equations 
of motion remain second order. 
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higher order nature of the field equations manifests itself and additional 
data, corresponding to a new dual operator, is required. 

A top down model arising from string theory must be consistent with 
unitarity. This is guaranteed if the curvature invariant is such that the 
resulting equations are actually second order. (Note that since one is treating 
the corrections perturbatively it is guaranteed that the shift to 7]/s is small 
and is consistent with unitarity, in contrast to the discussions of [8, 9] in 
which the coupling constant of the higher derivative term is allowed to be 
of order one.) 

As we discuss at the end of section 4 another case in which the higher 
derivative invariant is automatically consistent with unitarity is when the 
linearised field equation around AdS remains second order. This is a weaker 
condition than requiring that the equation of motion is always second order, 
but suffices to ensure that there is no non-unitary dual operator induced by 
adding the higher derivative term. Examples of such curvature invariants 
are those built out of the Weyl tensor of order three and higher. 

Finally it is interesting to note that reducing a curvature invariant of a 
given order from ten or eleven dimensions to four dimensions on a curved 
manifold gives rise to curvature invariants in the effective four dimensional 
action which is both of the same order and of a lower order, see for example 
(37). In the context of AdS solutions the reduction required is indeed al- 
ways on curved manifolds such as spheres. This implies in particular that a 
curvature invariant such as one quartic in the Riemann tensor never arises 
without an accompanying term quadratic in the Riemann tensor and a shift 
of the cosmological constant. Here we show that the term quadratic in the 
Riemann tensor gives rise to a new boundary condition for the linearised 
theory around AdS, and hence a dual operator in the CFT, which turns out 
to be non-unitary. When one combines all terms arising from the corrections 
at a given order in the upstairs theory, the resulting four dimensional theory 
must be unitarity and this may be achieved either by the linearised theory 
around AdS being second order or by the higher order terms conspiring to 
give a unitary dual operator. 

The plan of this paper is as follows. In section 2 we discuss in more 
detail the variational problem in higher derivative theories and show that it 
is well-posed with only boundary data for the metric when one treats higher 
derivative terms perturbatively about a leading order Einstein solution. In 
section 3 we first discuss what curvature invariants are expected to arise in 
the effective four-dimensional action from a top down perspective and then 
we explore the effects of various curvature invariants on four dimensional 
planar AdS black holes. Our goal is to find an analogue of the Gauss-Bonnet 
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corrected black hole in five and higher dimensions, i.e. a representative 
corrected AdS^ planar black hole, and we find that the solution in the Weyl 
corrected theory is the closest analogue. In section 4 we look in detail at the 
spectrum in theories with curvature squared corrections, demonstrating that 
there are indeed new dual operators associated with the higher derivative 
terms and these are non-unitary. Noting that the spectrum in the Weyl 
cubed theory is unchanged again this seems to be the simplest case of a 
representative correction. In section 5 we conclude. 

2 The variational problem in higher derivative the- 
ories 

In general one cannot define an analogue of the Gibbons-Hawking- York term 
[28, 29] such that the variational problem is well-defined with only a Dirichlet 
condition on the metric - one must impose additional conditions. This obser- 
vation explains a long standing problem in the literature: for generic higher 
derivative corrections the analogue of the Gibbons-Hawking- York (GHY) 
term has never been found. 

There is considerable literature discussing the variational problem in 
higher derivative theories. In the context of corrections arising in string the- 
ory, boundary terms were discussed in [31] where the analogue of the GHY 
term was found for Gauss-Bonnet. This is a very special case, however, as 
the field equations are second order. For corrections involving powers of the 
Ricci scalar, the variational problem was discussed in [32]. The generic issues 
in setting up a variational problem for higher derivative gravity given only 
a boundary condition on the metric were highlighted in [33]: the boundary 
terms which arise in varying the bulk action cannot in general be integrated 
to give an analogue of the GHY term. 

Here we argue that the problem in finding a GHY term results from the 
fact that in general such a term cannot exist: one must specify additional 
data together with the metric. In special cases an analogue of the GHY term 
was found, for example, for Lovelock theories. However, Lovelock theories 
are themselves special in that the equations of motion are actually second 
order and this fact explains why a GHY term could be found. 

A useful approach to dealing with higher derivative theories is the auxil- 
iary field method and the variational problem in such a context was discussed 
in [34]. In this approach the higher order equations are reduced to coupled 
second order equations for the metric and the auxiliary fields, and one spec- 
ifies boundary data for both the metric and for the auxiliary field. In the 
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context of perturbatively evaluating higher derivative corrections on leading 
order Einstein solutions, the boundary condition for the auxiliary field does 
not involve new data, but rather can be built out of the boundary data for 
the metric. When one looks at the spectrum, however, one sees that there 
is indeed generically new data required for the auxiliary field. These points 
will be illustrated further when we use the auxiliary field method to discuss 
the spectrum in section 4. 

Before moving on to consider specific models for higher derivative cor- 
rections in four dimensions, let us discuss the issue with the variational 
problem. We consider a general action in (d + 1) dimensions 

/ = / d d+1 x^C, (5) 
J M 

where the Lagrangian C depends only on the metric and the Riemann tensor. 
The variation of the action with respect to the metric gives 

SI = J M d d+l x^g-(^g^C + C^Sg^ (6) 

+ / d d+l x^~g U^R^ p x Sg aX + 2V p V M £^<w) 

+2 / azrc^vscr + ■■■ 

JdM 



Here dM is the boundary of the manifold M and we define 

while the ellipses denote boundary terms which vanish with a Dirichlet 
boundary condition on the metric, 5g = 0. 
In the case of Einstein gravity 

C = ^(R-2A) (8) 

and thus the boundary term in the variation is 

j d^(g^V^g ua -g^V u 5g^)- (9) 
JdM 

As is well-known, one can set up a well-defined variational problem by notic- 
ing that 

5 (-\ { dx^K] = ~ { dV{g™V»5g va + ■■■) (10) 

V K z JdM J 2K Z JdM 
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where the ellipses again denote terms which depend only the on restriction 
of the metric variation to the boundary (and which hence vanish given the 
boundary condition). Thus the addition of this term, the Gibbons-Hawking- 
York term, to the action gives a well-defined variational problem in which 
the metric on the boundary is held fixed. 

For generic Lagrangians involving higher powers of the curvature, the 
boundary terms involving metric derivatives cannot be canceled by those in 
the variation of a boundary term. To illustrate this it is useful to look at a 
specific example, 

C = ^(R-2A-aR pupa R^)- (11) 

When a = this reduces to the Einstein theory. When a is small the higher 
derivative term can be treated perturbatively. It is thus useful to express 
the equations of motion in the form 



Rjxv = T pu = Ag^ + atfiv, (12) 

4 

(d-1)" "^'"^ 1 ( d ~l) 
-2R W(rX Rf aX - 4V p V a R p ^ a , 



tfiu — /, 1 >-^' > ^ R parr\9 pv ~\~ -.-^p^/jR^ 9pv 



with T being the effective (trace adjusted) stress energy tensor. In later sec- 
tions we will be interested in four dimensional models, in which the Riemann 
squared term can be rewritten in terms of the Ricci tensor and the Ricci 
scalar, but in this section we will work in general dimension. The reason for 
consider terms involving the Riemann tensor is that such terms will always 
arise from top down models, and (unlike Ricci scalar and Ricci terms) they 
cannot be removed by field redefinitions. We can use the Bianchi identities 
to simplify the stress tensor as 

^■jxv = ^ 1)"^^ ^ Rp<JTr\9pv ~ 2R li p (J \Rj' (13) 

+^-^0^ - 4DiV + 4V^V MJ R, P , 



where □ = V p V p . 

A perturbative treatment of the field equations means that one looks for 
a solution such that 

9pv = 9{o)nv + ot9(i)tw H ( 14 ) 
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where is Einstein with cosmological constant A and satisfies 

(C R -A)g (1) ^ = j^^R paTn {g(Q))R P arr,{g(Q))g(Q)^ (15) 

-ZRppaX (9(o) ) R,y f " 7X (9(0)), 

where Cr is the linearized Ricci operator and terms on the right hand side 
are evaluated on the metric g/ \ using the connection of that metric. Note 
that the terms involve derivatives of the Ricci tensor do not contribute since 
the covariant derivative of the Einstein metric gv ) is zero. As emphasised 
earlier, this equation is a second order inhomogeneous equation for g^ and 
therefore it does not require any new boundary condition. Note that we 
regard here the boundary conditions for the metric as being given as a power 
series in a; i.e. the homogenous part of the solution g^ is determined by 
this data. 

Let us now turn to the question of the variational problem for such a 
theory. The new (relative to Einstein gravity) boundary term that arises in 
varying the action is then 

In r 

— d^R lu/pff VS S r + --- (16) 

where we again suppress terms which vanish for the boundary condition 
5g = 0. This term can be manipulated using the Gauss-Codazzi relations 
as follows. The metric on M can be decomposed as 

ds 2 = (N 2 + N^N^dr 2 + 2N^dx fl dr + i^daPda? (17) 

in terms of hypersurfaces S r of constant r with the unit normal to each 
hypersurface being given by n M . As the notation suggests, we are most 
interested in the case where the finite boundary is at spatial infinity, so r 
is indeed a radial coordinate 4 . Defining the radial flow vector r M such that 
r HQ r = ]^ components of tangent and normal to the hypersurfaces 
define the shift and the lapse Nn^ 1 respectively. The extrinsic curvature 
K^ v of the hypersurface is given by 

K^v = -Cnlnv, (18) 

where C is the Lie derivative. The Riemann tensor of the (d+1) dimensional 
manifold can now be expressed entirely in terms of the intrinsic curvature 

4 Such a foliation would also be appropriate near timelike infinity, in which case r would 
be a time coordinate and g rr < in Lorentzian signature. 
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and extrinsic curvature of S r via the Gauss-Codazzi relations 

^7u7p7 5 aRa(3yS = Rfiupa + K pa K up - K w K va ; (19) 
l> a R P a = DpKZ-DvKfr 

where D p is the covariant derivative of the metric 7 and R denotes the 
curvature of this metric. A useful manipulation of these equations gives the 
following identities 

K 2 -K^Kp V = R + ZG^nW- (20) 
£ n Kp U + KK tlu -2K^K pu = - ^ R p>7 , 

with Gp V the (bulk) Einstein tensor. One can simplify these expressions by 
fixing the gauge freedom such that N = 1 and N^ - = 0. In this gauge 

ds 2 = dr 2 + -f ij dx i dx j ; (21) 

Moreover the Gauss-Codazzi relations which we will need can be written in 
terms of the trace adjusted stress energy tensor as 

Rrirj = -d r Kij + KfK k f (22) 

K 2 — KijK^ = R + frr-flfij- 

drKij -2K\K lj + KK ij = Rij-fij. 

Returning to (16) the terms in the variation which do not vanish given a 
Dirichlet condition on the metric, $7 = 0, in this gauge take the form 



4/ d d x^R rirj 6^ = 4/ d d x^(-d r K il+ K?K kj )d r 5 7 lj (23) 

H z JdM « JdM 

= ^ f d d x^(KK tJ -K^Kkj-Rij + f^drS^, 
K JdM 

where in the last equality the bulk equation of motion in Gauss-Codazzi 
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form has been used. Next we note that 

5{KK l i Kij) = ^K i iK ijl kl d r 8 lkl + KK i id r 8 lij + --- (24) 

5(K 3 ) = Z -K 2 ^d r 6 lij + ---- 

S(RK) = ^Brfidrfryij + • • • ; 

5(K tJ K^ k Kl) = l K ^K)d r 5 llJ + ---- 

5(fr'Kij) = ^d^ij + ■■■; 

where ellipses denote terms which do not depend on the normal derivative 
of the metric derivation. Then 

^ f d d x^{KK l3 - K k K kj - Rn + I), '•' (25) 

K JdM 

= 3| f d d x^8(KK^K t] - \k z + RK- ^KfK^K^ - 2R ij K i >) 

K JdM & "J 

+ ^ / d d x^(^(ft - T rr ) + T^dr&yij +■■■ 
K JdM 

The terms in the second line are written in terms of quantities intrinsic to 
the boundary and define an analogue of the GHY term but the remaining 
terms left over in the last line cannot, in general, be expressed in terms of 
such quantities. 

Suppose however that one works perturbatively in a, evaluating the cor- 
rections as a perturbative series on the leading order metric, so that 

9w = 9$ + ctgty + ■■■ (26) 

Working to order a in the action one needs to evaluate the terms involving 
T only to zeroth order in a, i.e. 

(T fc fc - f rr ) + T*J ) -> dkfi (27) 

so that 

/ d d x^^(~f ij (ft -f rr ) + f ij )d r 5-iij -»• / d d x^^2dK5K + ■■■ (28) 

JdM JdM 

That is, applying the field equations perturbatively, the problematic term 
can indeed be reexpressed in terms of quantities which are intrinsic to the 
boundary. 
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Putting the terms together, we see that working up to order a the bound- 
ary term needed to set up a well-defined Dirichlet variational problem at a 
finite radial boundary is 



where implicitly in the first line one needs the metric to order a whilst in 
the second line one needs the metric only to zeroth order in a. 

It would be interesting to extend this proof to show that the variational 
problem is well-defined to arbitrary order. To do this one would need to 
argue that the problematic term 



can always be expressed as the variation of a term intrinsic to the boundary, 
when the bulk equations of motion are used iteratively. Such an all orders 
proof could be developed using similar inductive techniques to [35, 36]. 

3 Gravity models 

In this section we will consider higher derivative corrections to Einstein 
gravity with a negative cosmological constant in four bulk dimensions. Be- 
fore we describe the features of various models, let us comment on top 
down derivations of the effective action. One might think that it would be 
straightforward to work out the leading order corrections to the action from 
the reduction of ten or eleven dimensional actions, i.e. one could exploit 
our knowledge of the M theory action or the type II string actions. Here we 
point out that there are many subtleties in implementing such as strategy 
and our knowledge of these actions is not currently adequate to derive the 
corrections to AdS gravity actions in lower dimensions. 

To illustrate this point let us consider the best understood top down 
possibility to obtain AdS±, the reduction of M theory on a seven dimensional 
Sasaki-Einstein SEj to four dimensions. At the level of supergravity, it is 
always consistent to retain just the four-dimensional graviton in the lower 




K' 



% [ d d x^w\rt - f rr ) + 

' JdM 



(30) 
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dimensional theory, i.e. the eleven dimensional equations are solved by 
eleven-dimensional fields such that 



ds 2 u = ds\{E A ) + ds%(SE 7 ); (31) 
Fa = Ve 4 , 

where E4 is any Einstein manifold with negative cosmological constant, 
r7(Ki) is the volume form of this manifold and the metric reduction is di- 
agonal over the SEj. The effective four dimensional action is written only 
in terms of the metric on £4, g„ v . Note however that not only the eleven- 
dimensional metric g mn but also the four form F^ in eleven dimensions are 
non trivial, and the Riemann tensor of the SE7 is also non trivial since the 
manifold has positive curvature. 

Let us consider what this implies for the higher derivative corrections to 
the effective four-dimensional action. Since the four form is non-trivial at 
leading order, to compute the higher derivative corrections to the leading 
order solution, one would need to know higher derivative corrections to the 
eleven-dimensional action involving not just curvatures but also the four 
form. Building on [37] , [38] , leading corrections involving the latter in eleven 
dimensions were worked out in [39]; they have the structure 

I = J d u x^da (t 8 t 8 R 4 + ^eni? 4 ) (32) 

+ J d u x^g~b (t 8 t 8 R 4 - ^e u R 4 - ^e u t 8 AR 4 + [R 3 F 2 ] + [R 2 (DF) 2 ]^j , 

where a and b are coefficients. It is known by comparison with II A string 
calculations that 

6 = 2^2^47 Y' (33) 

with 2k\i = (27r) 8 /p. Here en is the eleven-dimensional epsilon, t 8 consists 
of 4 Kronecker deltas and t 8 t 8 R 4 denotes a specific product of such such 
t 8 tensors and four Riemann tensors; the explicit expressions will not be 
needed here. A is the three form of which F is the four form field strength. 
The tensor structure of the terms denoted [R?F 2 ] and [R 2 (DF) 2 ] is also not 
important here; all we need is this schematic form, in which D denotes the 
covariant derivative. 

One might think that the knowledge of such terms would suffice to com- 
pute the leading corrections to the eleven-dimensional solution of interest, 
(31), and that the these corrections could be rewritten in terms of a cor- 
rected equation for the four-dimensional metric g^ u , and hence in terms of a 
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corrected four dimensional action. Apart from the complexity of the actual 
calculation, there would be a number of subtleties in actually carrying this 
out. 

First of all, one cannot assume a priori that the higher order terms do 
not induce additional four-dimensional fields, as well as the metric, although 
it seems reasonable that in some cases they do not. For example, consider a 
four-dimensional massless scalar field 4> which corresponds to a modulus of 
the dual conformal field theory. In principle, even though this field is con- 
stant at leading order, it could be sourced by a higher derivative correction, 
i.e. one could have an equation such as 

U<j) ~ R n (34) 

where R n denotes schematically a scalar curvature invariant of order n. The 
latter must be zero when evaluated on AdS itself, as one does not expect 
the conformal invariance to be broken, but this argument could not exclude 
invariants of the Weyl tensor occurring. 

Even if could argue that a four-dimensional action involving only g„v 
exists, there is a second obstacle in actually computing such an action To 
illustrate this point, consider just one of the tensor structures occurring in 
the i? 4 invariant _ 

J^{RmnpqR Pq ) ■ (35) 
lp 

Evaluated on the lowest order metric this picks up contributions 

^ ((iW^"^) 2 + (R,u P aR^ Pa )(Ra h c d R abCd ) + (Rabc d R abCd ?) , (36) 
lp 

where R a bcd is the Riemann curvature of the Sasaki Einstein. For any given 
Sasaki Einstein this would then result in a term of the form 

/ ~ / d^x^-giiR^R^f + ^{R^R^) + bo), (37) 

in the four-dimensional action with Vs e 7 the volume of the Sasaki-Einstein 
and (60,62) computable (dimensionful) parameters. That is, a quartic in- 
variant in eleven dimensions can lead to quadratic and constant terms in 
four dimensions, with the latter shifting the cosmological constant. 

Since the Sasaki-Einstein has a curvature radius of the same order as 
the four dimensional manifold, none of these terms is subleading. Let L 
be the scale of the curvature radius for both; then each of the three terms 
in the action is of order L 3 /lp, using the fact that Riemann squared is of 
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order 1/L 4 . Note that the Einstein term in the action would be of order 
L 9 /lp. Terms arising from the reduction of higher order invariants in eleven 
dimensions would be subleading in a power series in L/l p . 

Recall that in the case of AdSi x S 7 the radius L scales according to 
L 6 ~ Nip where N is the rank of the dual gauge group. Therefore the 
Einstein term gives the well-known scaling of N 3 / 2 [40] whilst the terms 
given above scale as N 1 / 2 and thus are suppressed by a factor of 1/N relative 
to the leading order terms. We will use this scaling later when discussing the 
spectrum. Similarly for the case of ABJM [19] where the eleven-dimensional 
geometry is AdS^ x S 7 /Z^ the curvature radius scales according to L 6 ~ 
(kN)lp where N is the rank of the dual gauge group. Recalling that the 
volume of the compact space scales as 1/k this gives a scaling of k^ 2 N 3 / 2 
for the leading Einstein term. One can rewrite this as k 2 X 3 ^ 2 where the 't 
Hooft coupling is A = N/k, and this scaling was reproduced from the ABJM 
theory in [20]. In the ABJM case the term give above would contribute at 
order A 1 / 2 , i.e. it differs by a factor of l/(fe 2 A) from the leading term. (Note 
that validity of the eleven-dimensional description requires N k 5 .) 

In conclusion, identifying the leading order corrections in four dimensions 
is very subtle. The leading order correction in four dimensions indeed derives 
from the leading order correction in eleven dimensions, but terms involving 
higher curvature invariants in four dimensions can actually contribute at 
the same order! Similarly terms in the higher dimensional action involving 
R 3 F 2 and so on can give rise to corrections in four dimensions involving R 3 . 

Note that the term picked out above (35) shifts the cosmological constant 
and is non zero even when evaluated on AdS itself. This would mean, 
in particular, that it would be expected to adjust the value of the free 
energy for the dual theory (at zero temperature) evaluated on an S 3 . If one 
can argue that there is no such renormalisation, then the four-dimensional 
contributions from such a term must cancel those arising from the reduction 
of other eleven-dimensional terms. A series of corrections expressed in terms 
of the Weyl tensor, which vanishes on a maximally symmetric space, would 
not induce such a change in the free energy. 

From the string theory perspective one might think that one should in 
any case start from a higher dimensional with curvature corrections involving 
only the Weyl tensor, since corrections involving Ricci and Ricci scalar can 
always be absorbed into field redefinitions. Here we will look nonetheless 
look at terms such as (35) as well as Weyl terms. Firstly it is is interesting 
to look at the effects of different curvature invariant structures but secondly 
the usual field redefinition argument refers to the bulk field equations but 
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does not take into account boundary conditions and onshell thermodynamic 
quantities. We will see below that it is possible to have terms which do 
not contribute to the field equations perturbed around a given leading order 
solutions but which nonetheless change the action and change the spectrum. 
In particular, curvature squared corrections in four dimensions do not change 
the metric, so in the past they would have been viewed as trivial, but here 
we show that they still introduce additional (non-unitary) operators into the 
dual CFT spectrum. 

In what follows, we will pursue a bottom up perspective, in which we 
consider case by case the effects of various higher derivative terms in four 
dimensions. In other words, we discuss the effects of adding particular scalar 
curvature invariants to the four dimensional action. We will then return to 
the issue of which scalar invariants are expected to arise in top down models. 

3.1 Curvature squared corrections 

Motivated by requirements of renormalizability of gravity, curvature-squared 
modifications to Einstein's theory were first discussed in [41, 42] and they 
have been extensively explored in the literature. The most general action 
involves curvature squared terms can be written as 

1=^ J d A x ^(R-2A + aBTR^ + (3R 2 + T^^^W) (38) 
However, it is well known that the Gauss-Bonnet invariant, 

E 4 = R» upa R^ pa - IBTR,,,, + R 2 (39) 

does not contribute to the equations of motion in four dimensions but yields 
only a surface term. Hence for analysing the field equations we can eliminate 
the Riemann squared term in the action above, making the most general 
action we need to consider, modulo the E4 term, simply 

/ = ^2 / dA% V^9(R - 2A + aRTR^ + I3R 2 ), (40) 

(where implicitly the coefficients (a, j5) have been shifted relative to the 
above.) The equations of motion following from this action are 

+ E^ = (41) 

where 

Guv = Rnv - ^Rg^v + Ag^u (42) 
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and 

= 2a{R w RP-^RP°R pa g^) + 2pR{R^-^R gixv ) (43) 
+a{UR ilv + ^DRg^ - 2V p V^R p u) ) + 2p( 9lxl/ DR - V P V U R) 

By analyzing the above equations of motion it immediately follows that all 
solutions of the a = (3 = theory are also solutions of the full theory as E^ 
is zero for any Einstein spacetime. In particular, AdS-Schwarzschild black 
holes 

are solutions of the higher curvature theory. Here k = and e = cor- 
responds to the case in which the horizon is flat, with k = 1 and e = 1 
corresponding to the case in which the horizon is a two-sphere. Note how- 
ever that the thermodynamic properties are modified in the deformed theory 
and depend explicitly on the deformation parameters. 

It is straightforward to derive the thermodynamic properties in the de- 
formed theory, exploiting the fact that the metric remains Einstein. (For 
earlier discussions of thermodynamics in bulk dimensions higher than four 
see [43].) The free energy of the black holes can be obtained by consid- 
ering the onshell value of the action. In order for the variational problem 
to be well-defined, the action must be supplemented by boundary terms. 
For the Einstein part of the action the appropriate boundary term is the 
Gibbons-Hawking- York term discussed earlier 

Ighy = ~\ / d 3 xK^. (45) 

where K denotes the second fundamental form and 7 is the boundary met- 
ric. This term is not however sufficient to ensure a well-defined variational 
problem: in varying the bulk action the following boundary terms arise, 
analogously to those given in (16) 

\ [ d^ (2aR ua + 2(3Rg ua V fl 5g° u + • • • ) (46) 
K z JdM 

where we again suppress terms which vanish for the boundary condition 5g = 
0. Using the equation of motion R^ v = kg^ v we note that the variational 
problem will be well posed if we add the following boundary terms 

1 = -— J d 3 x^^K(2aA + 8/3 A). (47) 
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In this case the fact that the solution remains Einstein implies that this 
term is sufficient to evaluate the onshell action, to all perturbative orders in 
a and f3. (It does not however suffice for discussing fluctuations around the 
Einstein solution, as we will discuss in section four.) 
Evaluating the complete onshell action gives 

I = -L 2A(i + 2aA + 8/3A) J d 4 x^/^g (48) 

-\ (1 + 2aA + 8/3 A) f d 3 x^jK. 

Relative to the case of a = /3 = 0, there is just an overall prefactor, which 
means that we can immediately read off from [40, 44, 45] the required holo- 
graphic counterterms as 

la = -2^2 (1 + 2«A + 8/3A) J d 3 x^ (± + Zi?( 7 )) , (49) 
where Z 2 = |3/A|. The asymptotic expansion of the metric g is [40, 45] 

ds 2 = I 2 + \g ij (x,p)dx i dxA ; (50) 



ftj (x, p) = £/ (0)ij (x) + p 5(2)^- (x) + p g {3)ij + ■■■ , 

ff(2)ij = -Rij{g(p)) H ^ — 5( )ij, 

with 5(3) being traceless and divergenceless but otherwise undetermined by 
the asymptotic analysis. The renormalized stress energy tensor obtained by 
varying the action with respect to <7( ) is then shifted by an overall prefactor 
relative to [45] 

< T u> = A (1 + 2«A + 8pA)g m . (51) 

We can now immediately evaluate thermodynamic quantities for the black 
hole solutions (44); the free energies and masses are clearly shifted relative 
to those in Einstein gravity by a proportionality factor: 

-frF = /f nshell = ^|I(l + 2«A + 8/3A)m (52) 

M = ( d 2 xyf^(T m ) = + 2aA + 8/3A)m (53) 

with 0t the inverse temperature (not to be confused with the coupling con- 
stant (3) and the temperature being 

2 I A I m \ 

-Yr h + ^\, (54) 
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and rh is the horizon position. I E denotes the Euclidean action, which in 
this static case is straightforwardly computed by analytic continuation of 
the time. Note that under such a continuation il — > —I E . One can also 
work out the black hole entropy using Wald's method [46]. Define 

qv, = _2C^ pa V p l a + --- (55) 

where Q^ v is antisymmetric and the terms denoted by ellipses vanish for 
stationary horizons. Here 

C^ pa = JL . (56) 
For a stationary horizon the black hole entropy is then given by 

S = ^j n Q? v d^, (57) 

with T being the horizon temperature, % denoting the horizon and l K being 
the horizon normal. 
For Einstein gravity 

L = ^^-g(R-2K) (58) 

where k 2 = 8ttG and G is the Newton constant. Hence 

SL y/=g 



^ Pg ua_ g ^ g u P) (5Q) 



And so 



and 



f n V=g-(WdX^) = ^ (61) 



i r , , A h 

with Ah the horizon area, using l u V v l p ' = where Kh = 2irT is the 

surface gravity of the horizon. 

For the curvature squared corrections, using 

xr>2 

' - = fl(f«r-ff); (62) 



^Rfiupcr 



5(R TT >R Tr >) = (RW a av - R^a vp ) 



20 



the Wald entropy becomes 



S= (l + 2aA + 8/3A)-^. (63) 

Putting these results together one finds that the thermodynamic relations 
F = M -TS; dM = TdS, (64) 

are indeed satisfied. 

To summarise, the metric is uncorrected but the thermodynamic proper- 
ties of the black holes are adjusted: the entropy, the temperature, the mass 
and the free energy are all changed, albeit by just an overall factor. It is also 
interesting to note that the action evaluated on AdS± with S 3 boundary is 
also changed. In the latter case the relevant bulk metric is 

ds 2 = — [dp 2 + sinh 2 p dttf) , (65) 

where < p < oo. Using the renormalised action given above, one can 
compute the onshell Euclidean action to be 

127T 2 

Onshell = (1 + 2A« + 8A/3). (66) 

This is therefore corrected by the curvature squared terms except when 
a = —4:{3, which corresponds to the case in which the correction is Riemann 
squared minus E±. Given a holographic dual in which one can compute the 
free energy on S 3 by localisation techniques, the answer will give a criterion 
restricting the terms which can arise in the effective four-dimensional action. 
In particular, the case of ABJM theory, for which the exact expression for the 
planar free energy was obtained in [20], will be explored in detail elsewhere. 

At this point it would seem as if the addition of such terms to the action 
is rather trivial because the thermodynamic quantities are shifted by an 
overall factor, which could be reabsorbed into the cosmological constant. 
However, we will discuss in section 4, these terms are highly non-trivial 
when one looks at the spectrum of the theory. To find the spectrum of the 
dual CFT linearize the above field equations about the background solution 
AdS±. As we discuss in section 4, the bulk theory is found to describe a 
massless spin-2 graviton, a massive scalar and a massive spin-2 field. By 
tuning the coefficients so that a = — 3/3 one may eliminate the massive 
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scalar mode. One can also tune the remaining coefficient f3 to the so-called 
critical value [22] 

> = -u < 67 > 

where the massive spin two mode becomes logarithmic [47] . Noting that the 
AdS-Schwarzschild mass when a = —3(3 behaves as 

Tfl 

M = — (1 + 2/3A) (68) 

we see that in the critical theory the black hole solution has zero mass. One 
can show furthermore that the Wald entropy vanishes at the critical point. 

However we should emphasize that the critical value (67) can clearly 
never be achieved when the Planck length is small compared to the curvature 
radius of AdS. If one is viewing the higher curvature corrections as arising 
from a top down string model then (5 'A is necessarily much smaller than one. 
The critical theory does not therefore provide a good model for corrections 
to macroscopic AdS^ black holes. 

Let us now make a connection to conformal gravity. With the first 
parameter choice of a = —3/3 we may rewrite the higher curvature term in 
terms of the Weyl tensor: 

_I a(jR 2 _ 3R ^ Rfiu) = ± a (C^°C^ pa - E 4 ) (69) 
Hence the Lagrangian is equivalent to 

1 = i / * x ^ (R - 2A + \a(C^C pvp(T - £7 4 )) • (70) 

By taking the limit of a — > oo one recovers Weyl gravity, see related discus- 
sions in [48], but again this would not be reached as a small correction from 
an Einstein solution. 

The variation of the Weyl squared term with respect to the metric is lin- 
ear in the Weyl tensor. This means that this correction vanishes identically 
when evaluated on AdS^ since it has vanishing Weyl tensor. Therefore we 
can deduce from (66) that the onshell Euclidean action evaluated on AdS^ 
with S 3 boundary for 

I E = / d ^ x ^9 {R-2K- ^aE^j (71) 

is 

E 12vr 2 2qA 

-'onshell — pq^2 ( 3/" 
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In other words, the topological invariant does of course contribute to the 
action even though it does not affect the field equations. The renormalised 
E4 term captures the Euler invariant of the manifold with S* 3 conformal 
boundary. Tuning to the critical value (67) this action is zero. 

3.2 f(R) Gravity 

In our exploration of corrected black hole solutions we will now move on 
to consider an f(R) theory. The f(R) theory is obtained when we add a 
generic polynomial in the Ricci scalar R to the usual Einstein action, 

1 = ~h I d4x ^ ( R - 2A + f( R )) ( 73 ) 

where 

f(R) = ]T a n R n , (74) 

n>2 

with arbitrary coefficients a n . There us considerable interest in f(R) the- 
ories in the context of phenomenology and cosmology, see the reviews of 
[49, 50, 51], even though such corrections are not well motivated from top 
down considerations, since they can be removed by field redefinitions. 

Indeed it is well known that such a correction will not change the leading 
order black hole solution non-trivially, although it will change its thermo- 
dynamic properties. This follows from the equations of motion 

Guu, + = (75) 

with Gfiu defined in (42) and 

= J2 a n nRn ' 1 (R^-^Rg^u) (76) 

n>2 

+ £ a n n(^Di?"- 1 - V^ v R n - 1 ). 

n>2 

Consider an Einstein solution g^ v which satisfies 

Riiv = R^u = -^Rg^u, (77) 

Evaluated on such a solution both the second line together with the n = 2 
term in the first line of F^ v vanish and 

F^ = ]T a n (n - 2)(4\) n ~ 1 \g^ = -SAg^. (78) 

n>2 
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The field equations (75) are then satisfied provided that 



A = A + ^a„(2-n)(4A) n - 1 A; 

n>2 



(79) 



i.e. the higher derivative term acts to shift the effective cosmological con- 
stant. Treating the f(R) term as a small perturbation around the leading 
order solution by setting all coefficients a„ < 1, we may express 



A « A(l — a 3 (4A) 2 - 2a 4 (4A) 3 + 



(80) 



where in the non-linear terms we use the leading order behavior A ~ A. 
Therefore any Einstein solution remains an Einstein solution in the corrected 
theory, but with a shifted cosmological constant. 

One should again note that the corrected theory does admit non-Einstein 
solutions, but any solution which reduces to an Einstein solution in the 
leading order theory remains Einstein in the corrected theory. In other 
words, when one treats the higher order terms perturbatively one discards 
solutions which do not reduce to Einstein solutions on setting a n to zero. 
The higher derivative terms with recur, however, when one discusses the 
spectrum as one will obtain new propagating modes. 

As a warm up exercise for the non-trivial corrections discussed in the 
following sections it is useful to derive the corrections to static solutions as 
follows, using a similar method to that of [52] and also [53]. Let the metric 
be parameterized as 



dv 

ds 2 = -a(r)b 2 (r)dt 2 + ^— + r 2 (dx 2 + dy 2 ), 

a(r) 



(81) 



where we now focus on the case of flat horizons. Substituting this metric 
ansatz, the action reduces to 



1 = 2^ pTVxyI 



(82) 



with @t the periodicity in time; V xy the regulated volume of the (x — y) 
plane and 



/ 



dr 



A A n 

2 ^(-r 2 - a)+ B- 1)na «7^ 



-Abr 3 + 2ab'r 2 + a'br 2 + 2abr 



(83) 
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where we have used the fact that the Ricci scalar is given by 

b{r)r z 

with A(r) given by 

A(r) = 3a'b'r 2 + 2ab"r 2 + a"br 2 + 4a'br + Aab'r + lab. (85) 

The second line in (83) arises from partial integrations. Varying the bulk 
term in the action we find the following equations of motion for a and b: 

= 2rt '-i:r-(-ir o „A(^) (86) 

where 

5 f A n \ nA n - x , ,.. 9 . d ( nA n - 1 , ». A 

d 2 ( nA n - x „ 2 .\ 

The other equation of motion is 

5 / A n 



= -2Ar 2 - 2a - 2a'r - V l n -\-l) n a n ^- (- 

„ So V( 



! 56 V^r 2 )™- 1 
where 

5 / A™ \ (1 - nU n nA™" 1 , „ 9 

»((S3pi) = W^ + m^ (<1 ' +4o+2o) (89) 

d / nA™- 1 2 , , s \ d 2 ( nA n - x X 



d7 1 Wr~^ r a + 4ar) J + I (6^ (2ar 



Solving these equations of motion to linear order in the coupling constants 
ttj we expand as: 

a{r) = O( ) (r) + ^a n a (n) (r); (90) 

n 

= &(0)( r ) +H«n6(n)( r )- 

n 

To leading order, namely all a n = 0, the equations are solved by 

6 (0 )(r) = l; a {0) (r) = -Ur 2 -^ (91) 
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To linear order in the perturbations, the general solution to the equations 
of motion is 



6 (n) (r) = 0; (92) 
a ( „)(r) = (-f(n-2)(4A)«- 1 Ar 2 -^. 

The latter renormalizes the cosmological constant and in addition allows for 
a shift in the integration constant which parameterizes the black hole mass: 
the corrected metric is 

ds 2 = —r— — - dt 2 (--Xr 2 - - - V ^^) + r 2 dx ■ dx. 

(93) 

3.2.1 Black hole thermodynamics in f(R) theory 

For the f(R) term the analysis is of the variational problem is subtle: varying 
the bulk term gives rise to a boundary variation 

8I = \I d 3 xf'(R)5(K^), (94) 

where f'(R) = dnf(R). The appropriate boundary term for a four-dimensional 
f(R) theory was argued by Hawking and Luttrell [32] to be 

Ihl = ~ J d 3 xf'(R)K^ (95) 

However, in general this is not satisfactory since R is not intrinsic to the 
boundary, i.e. it is the scalar curvature of the bulk metric, rather than the 
boundary metric. In the case at hand however one can use the fact that 
the onshell Ricci scalar is constant to write this term in terms of quantities 
manifestly intrinsic to the boundary. Putting all terms together the complete 
action is 

1 = 2^ Jd 4 x^(R-2A + f(R)) (96) 
__L J d 3 x^K(l + f'(R)). 

To evaluate the free energy one needs to holographically renormalize this 
action. To linear order in the couplings for the higher derivative terms, 
however, one can immediately carry out this procedure using the known 
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results for asymptotically locally AdS Einstein manifolds. Recall that the 
deformed solution is Einstein, with a different cosmological constant, the 
metric can always be expressed as 

ds 1 = Pds 2 (97) 

with 

P = pr; V = ( 98 ) 

and A < 0. The asymptotic expansion of the metric g is known 

ds 2 = -At- H — T7(jij(x, p)dx t dx : > '; (99) 

<7ij (s, p) = 5(0)ij 0*0 + P 2 5(2)ij (^) + P 3 <?(3)ij + ■ ■ ■ , 

j? r- n i ^(g(Q)) - 
ff(2)ij = -MijK9(p)) H ^ — 5(o)ij> 

with <7( 3 ) being traceless and divergencless but otherwise undetermined by 
the asymptotic analysis. Given this form for the asymptotic expansion one 
can now compute the regulated action and hence the counterterms. In doing 
so one can use the fact that to linear order in the new couplings 

f(R) -> f(R)\ R=iA . (100) 

Setting A = — 3 so that the metric to leading order is normalized to unit 
curvature radius, the required counterterms are then 

let = "2^2 (l + E 2 ««(- 12 ) n_1 ) Jd 3 x^f(4 + R^)); (101) 

= ^(i + E(Y- 1)a " ( - 12)n_1) / d3x ^ (4) 

+ 2-2 (1 + E( Y - 1 )«-(- 12 ) n " 1 ) 7 <*V=^(1 + £(n - ^^(-l^- 1 )^), 



and the renormalized stress energy tensor obtained by varying the action 
with respect to is 

< T ^> = ^(1 + E 2 ««(-12) n - 1 )5(3) ir (102) 
One can then compute the mass of the black hole in (93) as 

M = jd 2 x(T 00 ) = ^|(l + ^2a n (-12) n - 1 )(m + E«nm n ), (103) 
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and evaluating the onshell action gives 



-frF = I E ee= /3 T ^|(1 + £ 2a n (-12)^ 1 )(m + £ a n m n ) (104) 

n n 

with F the free energy and the black hole temperature being 1//?t- 
In the /(i?) theory using the fact that 

r) R n 1 

nR n -\gWg™-g^g v P), (105) 



^RfMUpa 2 



the Wald entropy (57) is given by 



971-4, 

S={l + "£a n (U) n - 1 )—^. (106) 

n K 

Evaluating this one obtains 

S = ^(1 + J2 2a n (-12)»- 1 )(m + £ a„m n ) 2 / 3 . (107) 

K n n 

Finally the temperature of the black hole is given by 



T=— (m + ^onron) 1 / 3 . (108) 



Putting these results together one sees that the relation F = M — TS is 
satisfied together with the first law dM = TdS. Moreover, it is clear that 
by choosing the integration constants m n such that 

m n = -2(-12) n_1 m (109) 

the black hole in the f(R) theory has unchanged thermodynamic properties 
to leading order in the coupling constants a n . 



3.3 Einstein + C 3 

We now move on to consider the addition to the action of curvature invari- 
ants of degree three or higher. At this point it is useful to look at classi- 
fications of scalar curvature invariants in our dimensions. One such set of 
invariants are the Carminati-McLenaghan invariants, [54]. At degree three 
the possible invariants include both those built of lower degree invariants, 
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for example RC^^C^pa, and invariants built by contracting three tensors 
with each other. At degree three the latter gives the new invariants 



° D pi" :, p 

qpv q (~l 



ripv f^P® riTr\ 

Spa * Cppua 



fipv f-tpo (-ITT) 
* U po^ TT\ fMU 



(110) 



where Sp_ u is the traceless Ricci tensor, Cp Vpa is the Weyl tensor and *Cp Vpa 
denotes the dual of the Weyl tensor. Note that these comprise an over 
complete set of invariants for a planar static spacetime. Curvature invariants 
built from the Ricci scalar or Ricci tensor will behave qualitatively similarly 
to those at quadratic order, leaving the metric unchanged but shifting the 
action. Therefore in this section we will focus on the effect of the cube of the 
Weyl tensor on planar black hole solutions, which is qualitatively different. 
The action we consider is therefore 

1 ' ' -d(R-2A + aC fM rC p J' X C v r) (HI) 



2k 2 



d 4 x 



where Cp Vpa is the Weyl tensor. 

Since the general field equations are somewhat complicated in this case, 
the easiest way to obtain the corrections to the planar black holes is as 
follows. Evaluated on a static ansatz (81) the action reduces to and effective 
one dimensional action 



= dr 
Jr H 



A 2^ B 



-Abr 3 + 2ab'r 2 + a'br 2 + 2abr 



rh 



where 



B = 2ab- 2rbo! - 2rab' + 3r 2 a'b' + r 2 ba" + 2r 2 ab" 



(112) 



(113) 



Varying this action we get the equations of motion for a and b and these are 
solved perturbatively in a as in Sections 3.2 and 3.4. There is in this case a 
non-trivial correction to the planar black hole solution: 



8m 2 A 16m 3 



and 



m 



(114) 



(115) 



Here a\ and b\ are arbitrary integration constants. The former acts as a 
redefinition of the mass parameter m at order a and the latter changes the 
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norm of the time Killing vector at infinity at order a. We will discuss the 
interpretation of these integration constants further below but setting them 
to zero we obtain 

a i r ) = «(0)( r ) +ao(i)(r) 

A 9 m ( 8m 2 A 16m 3 \ 

= ~r '-T +a [-—-—)- <U6) 

H r ) = &(o)(>") +a6(i)(r) 

/ m 2 \ 
= 1 + « ^-6^6- J ■ 

Note that the AdS solution itself is uncorrected, as one would expect: the 
contribution to the field equations from the variation of the Weyl cubed 
term is given below in (122) and evaluated on a solution with vanishing 
Weyl tensor it is zero. 

3.3.1 Thermodynamics of corrected black hole solutions 

Let us now work out the thermodynamics of the corrected black hole so- 
lution. The horizon is given by th such that a(ru) = 0. Since a(r) = 
a (o)( r ) + aa (i)( r )) we fi n< i r H also to order a. Let 

fH = r H (o) + ar H( i) (117) 

where 

rU = ~; (118) 



H(0) ~ X' 

r H (i) = ^(0)1 -33 A 



2 



The temperature of this black hole solution is given by: 



T = v ' v ' (119 

47T 

|A 



: r H(o) (l 

We can also work out the black hole entropy (57) giving 



4vr— r ° A 27 



- \2 



T .m(i + (120) 

AG 27 v ' 
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Note that although both the temperature and the entropy are corrected at 
order a the combination TS is actually uncorrected at this order. Moreover, 
imposing the thermodynamic relation 

dM = TdS, (121) 

and using the fact that the C 3 term evaluated on pure AdS is zero, we 
can infer that the mass must also be unchanged at order a. (In varying 
the entropy note that both a and A are held fixed.) Using the relation 
F = M — TS we can also then infer that the onshell action must also be 
unchanged at order a. 

One can also argue that the free energy and mass are unchanged at order 
a by considering their direct evaluation. Let us consider first the onshell 
action. The first step is to ensure that the variational problem, including 
the additional C 3 term, is well-posed at finite radius. To investigate this 
we vary the bulk action (111) with respect to the metric. From the term at 
order a one obtains the following contribution to the bulk field equation 

puv 1 „u,v /~i Triri Xk/^i po" ar 1 lir\/~iv piXnoa , o/^ttiXk/^ per r>v 

y ~ 2 P a TT ) ~~ DL/ pcr U r,XK^ + 6Lj U Ak U vrri 

+ AC tvXk C Xk ^R T(T - 2C TvXK C x ^Rg Ta - C^C parri R^ 

+ QV a V T {C TvXK C x ^) + 2V T V p (C T ^ XK C XK pu g pv ) + (122) 

-a^c^c^gpn + c^c^gn - 2V a v T (c Tr > XK c x pa g PV g fMU ) 

+ 2V a V^C^ XK C x pa gpv ) + V^V^C^C^), 

where Q^ v = R pu — \Rg iJLU + Ag^. The variation results in the following 
boundary terms at order a involving derivatives of the metric 



K Jd 



dM 1 

+ n x C p ^C n r9^Mx P ~ n c C^ T C v rg^^g pu 
- npC^C^rg^g^u^x - C^C^ pa SIC ] , (123) 

with n the normal to the boundary. (There are additional boundary terms 
involving the metric which automatically vanish for a Dirichlet boundary 
condition.) As one would have anticipated, boundary terms involving the 
normal derivative of the metric arise in this variation and it would therefore 
seem as if one needs additional Gibbons-Hawking like terms in order for 
the variational problem to be well-posed. Moreover, working iteratively in 
a and then using the bulk field equations to simplify the boundary terms 
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looks a very non-trivial calculation in this case. However, it turns out that 
one only needs to use the fact that the leading order metric is Einstein and 
is asymptotically locally AdS: in Feffer man- Graham coordinates (50), the 
leading power in the Weyl tensor necessarily behaves as 

<W ~ 1 (124) 

One can use this behaviour to show that the boundary terms needed for the 
variational problem to be well-posed all go to zero as a positive power of p. 
For example, the term 

\ [ Sxy/^C^C^K (125) 

K JdM 

evaluated at p = e <C 1 behaves as e, and therefore does not contribute 
in the limit e — >■ 0. Therefore, although one could indeed use the explicit 
expansion of the onshell Weyl tensor to express the boundary terms in terms 
of quantities intrinsic to the boundary, the resulting boundary action cannot 
give a finite contribution to the onshell action. 
The onshell action is thus given by 

/onshell = \ I d 4 x^(A + aC^C^C^) (126) 

~ J d 3 x^ (k + 2/1 + l -R(jf) , 

where we have used the onshell relation R = 4 A + C 3 . The terms in the 
second line denote the Gibbons-Hawking term along with the counterterms. 
The latter suffice to remove the divergences at leading order, but do not 
in general suffice to remove additional divergences at order a. However, 
it again turns out that the Weyl correction falls off sufficiently fast at the 
boundary that there are no additional terms needed at order a. To show 
this one needs to use the fact that, in Feffer man- Graham coordinates, C 3 is 
of order p e or smaller and the correction to the metric at order a is of order 
p 6 or smaller. Looking at the terms in the onshell action, this means that 
the contributions at order a are of order e 3 or smaller, and thus vanish in 
the limit e — > 0. For example, the term 

J d 4 x^C 3 ~ jdp^-p*- e 3 . (127) 

Thus only the above terms are needed in computing the renormalised ac- 
tion. Note that this implies that, as expected, the onshell action for AdS 
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is uncorrected at order a, regardless of the choice of conformal class of the 
boundary metric. In particular, if one computes the free energy for the dual 
theory on an S 3 , it is not changed at order a. 

It is still non-trivial that the actual value of the free energy for the planar 
black hole is uncorrected, as the metric is corrected, the horizon position is 
shifted and the C 3 term in the action all give finite contributions at the 
horizon. Explicitly evaluating the onshell action using (116) together with 



12m 3 



c 3 = T < 128 ) 



one obtains 



F = -^ nsh eii = \ rdrr%(r)(K + a 1 ^) (129) 

K Jr h V 

where r c 3> 1 is the cutoff radius. Integrating the bulk term and looking at 
the horizon contribution one obtains 

V xv ,2 k o 4am 2 A 4am 3 . V xv 

-f -Ar 3 + 3 —) rh = -f 2m. 130 

k z 3 r 6 r° k z 

i.e. the terms of order a cancel! Looking at the contribution from the cutoff 
boundary, as already argued the terms of order a fall off too quickly to 
contribute and one is left with a contribution 

(13D 

with the total free energy being 

F = ~^m, (132) 

i.e. unchanged at order a. One can similarly argue why the mass M = 
V xy m/ n 2 is unchanged at this order: varying the renormalised onshell action 
with respect to the source for the stress energy tensor, all terms at order a 
are subleading in the radial expansion and do not contribute. 

To summarise: the C 3 term leads to a correction of the metric of the 
planar black hole. The temperature and entropy are both changed at order 
a but the mass and the free energy are unchanged. 
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At this point we return to the physical interpretation of the integra- 
tion constants in the corrected solution. The first integration constant a\ 
corresponds to a shift in the mass parameter, 

m— >m — aa±. (133) 

This shift will affect the entropy, temperature, mass and free energy. The 
second integration constant corresponds to a redefinition of the time coordi- 
nate and hence of the temperature. One can see this by looking at the form 
of the metric 

ds 2 = -(l + 2t*6i + ---)a(r)dt 2 + ■■■ (134) 
= -a{r)dt 2 + ■■■ 

i.e. by redefining the time coordinate one can absorb the integration con- 
stant 61. This in turn corresponds to a shift of the temperature by 

T -> T(l + a&i), (135) 

with the free energy and mass shifted by the same factor. 

By an appropriate choice of the integration constant a± = —mA 2 /9 one 
can make the entropy be uncorrected at order a. However, this value of a\ 
is such that the temperature, mass and free energy are corrected: 

ah? ah? ah? 

T->T(l-2__); M ->■ M(l — ); F — >■ F(l — ). (136) 

(These corrections are clearly consistent with the thermodynamic relation.) 
By fixing the integration constant b\ appropriately and redefining the time 
coordinate, one can undo these corrections at order a, leaving all thermo- 
dynamic quantities unchanged to order a 2 . However, such a redefinition is 
somewhat unnatural from the perspective of the holographic duality, as it 
implies that the time coordinate for the field theory is redefined at order a. 

Thus the thermodynamics at order a depends on which quantities one 
has chosen to hold fixed. In the context of supersymmetric black holes one 
fixes the mass (and charge), with the temperature necessarily being zero and 
the entropy being corrected by the higher derivative terms. In the context 
of finite temperature black holes, it would seem natural to fix the mass also, 
as we did above, with the temperature and the entropy being corrected. 
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3.4 Einstein + R 4 



In order to obtain an eleven-dimensional correction which cannot be ren- 
dered trivial by field redefinitions we need to add a curvature invariant in- 
volving the Riemann tensor, with the first non-trivial term arising at fourth 
order. As discussed earlier, the reduction of this term will result in terms 
quartic in the Riemann tensor in the effective four dimensional action. In 
this section we will work with one representative curvature invariant at this 
order, the same tensor structure considered earlier, but the generalization 
of the analysis to other tensor structures would be straightforward. The 
reason for considering this particular term is because, we discussed earlier, 
such a term would accompany any Riemann squared term occurring in the 
effective lower dimensional action. 
The action we consider is therefore 

/ = ^2 / d4x ^9 (R-2A + a{R parX RP° TX ) 2 ) (137) 



The resulting field equation is 

1 

2 



G» v = l(R P arxR paTX ) 2 9^ ~ AR parX R^ x R^' 5 R^ 5 (138) 



+ 8V p V a (Rr XlS R T ^ s R^n- 

Note also that there is a new boundary term involving derivatives of the 
metric variation obtained when varying the term at order a 



K 2 



/ d 3 x « R p arxR p(TTX R aM& n 7 Vp6g aS . (139) 

JdM 



We will discuss this term in the context of the variational problem below. 

Now let us turn to the effect of such a correction on the planar black 
hole metric. Evaluated on the static ansatz (81) the action reduces to 



L 



oo 

dr 



D 2 

2rb'(a - r 2 ) + a- 



(140) 



(br 2 ) 

where to simplify formulae in this section we have imposed A = — 3 and 

D = 8a 2 b' 2 r 2 + Saa'b'br 2 + 4a' 2 b 2 r 2 + Aa 2 b 2 + 9a ,2 b' 2 r 4 + 12aa'b'b"r 4 
+6a'a"b'br 4 + 4aW + 4aa"b"br 4 + a" 2 b 2 r 4 (141) 

Varying this action we obtain equations of motion for a and b consisting 
of the piece coming from the Einstein action and a piece proportional to a 
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coming from the correction. Once again we solve these perturbatively in a 
as in Section 3.2. The planar black hole solution is indeed corrected: 

a ( r ) = a (o)( r ) + ( 142 ) 

«i nr . 2 \ 672m 2 1200m 3 536m 4 \ 
96r — 




r 4 r 7 r 10 



and 



b(r) = b (0) (r)+ab {1) (r) (143) 
/ 336m 2 224m 3 \ 



1 + a 6i + = h 



i y6 ^9 / 

Again there are two arbitrary integration constants, corresponding to re- 
defining the mass parameter and the time coordinate at order a. These 
constants will be set to zero and their effect will be discussed further below. 
Note however that the AdS metric itself is corrected by the quartic Riemann 
term, since the latter does not evaluate to zero, unlike the previous Weyl 
example. 

3.4.1 Thermodynamics of corrected black hole solutions 

Let us first calculate the corrected horizon position and temperature of this 
black hole solution. The horizon is located at 

nr = -g-(l + a-g-), (144) 
and the temperature of this black hole solution is given by: 

T = -^-m 1 / 3 (3 - 208a) (145) 

We also work out the Wald entropy using (57) giving 
A 4a 



S =40 



--f / (R"-"-R in , p ,)R"""^d i x, (146) 

K J H 

where we use ingoing coordinates for the horizon, namely 

ds 2 = -a{r)b{r) 2 dv 2 + 2b{r)drdv + r 2 (dx 2 + dy 2 ), (147) 

and the integral is over the spatial part of the horizon. The integrand in the 
second term however vanishes when evaluated on the leading order solution 
since 



1 1m 

R vrvr = i // = 2 _ till ( 148 ) 

2 r 3 
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which vanishes at the horizon. (The Riemann squared term is non-zero at 
the horizon.) Hence, in this case, the black hole entropy is given only in 
terms of the corrected area of the horizon 

S=^ m ^{l + Jf). (149) 
Combining the entropy and the temperature we notice that the combination 

TS =^m + 0(a 2 ) (150) 

is again unchanged at order a. As in the previous section we can now argue 
that for the relation 

dM = TdS (151) 

to hold the correction to dM at order a must also vanish. 

This argument on its own however does not exclude there being a term in 
the mass (and free energy) at order a which is independent of the parameter 
m: recall that, unlike the Weyl example, the AdS metric itself is corrected 
and the R A term in the action evaluated on AdS is non-zero. In other words, 
the holographically renormalized higher derivative action evaluated on AdS 
could be non- vanishing. Computation of the renormalised mass and action is 
somewhat involved as it requires analysing the corrections to asymptotically 
locally AdS solutions, isolating the divergences, computing the counterterms 
and so on. Fortunately there is a short cut: when the dual field theory is 
supersymmetric, the mass of the m = solution is necessarily zero, as is 
the free energy, and therefore there cannot be any contributions to the mass 
and free energy at order a which are independent of m. (Note that the 
free energy of the dual theory on a curved space would in general indeed be 
expected to be corrected at order a.) 

Thus, in summary, as for the C 3 case, the temperature and the entropy 
are corrected whilst the mass and free energy are not. By choosing the 
integration constants a\ and b\ one can adjust which thermodynamic quan- 
tities are corrected at order a, but the most natural physical choice from 
holographic considerations is indeed that where the mass is fixed and the 
entropy is corrected. 

3.5 Corrections arising in string theory 

In this section we have explored the effect of various curvature invariants 
added to the four dimensional action. We have shown that the thermody- 
namic properties of black hole solutions are in general corrected even when 
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the metric is not corrected. From a top down perspective it would be compli- 
cated to determine which curvature invariants arise in the four dimensional 
action, with a given higher dimensional invariant contributing to invariants 
of different derivative order in four dimensions. 

From the dual holographic perspective, one can try to restrict the in- 
variants which arise in four dimensions using the free energy on an S 3 . This 
would not restrict at all Weyl invariants which do not contribute to the 
free energy. One would also anticipate that other specific combinations of 
invariants involving Riemann, Ricci and Ricci scalar can be made in which 
the correction to the free energy also vanishes. 

In the following section we will turn to another criterion for higher deriva- 
tive corrections: the spectrum of fluctuations and the corresponding dual 
operators. In general, imposing that such corrections lead to CFT opera- 
tors which are unitary and have positive norm, rules out many curvature 
invariants. 

4 Spectrum of curvature squared theories 
4.1 Linearized equations of motion 

In this section we discuss the spectra of the higher derivative theories. For 
the sake of brevity, we will mostly focus on the case of curvature squared 
corrections, but the analysis for other higher-derivative theories would be 
similar and will be discussed at the end. 

We consider again the action (40), whose equations of motion are given 
in (41) to (43). Since our interest here is in the context of holography, we 
consider the spectrum of excitations around AdS^, and when we need an 
explicit form for the metric we will work in the Poincare patch in which 

ds 2 = ^- + \(-dt 2 + dx 2 + dy 2 ). (152) 
p z p z 

We denote by x l the coordinates on the three dimensional slices of constant 
p. We vary the metric as g^ u ->■ g^ u + Sg^ = g^ + h^, where g^ u is the 
AdSi background metric. It was shown in [22] that the linear variations of 
the various tensors appearing in the equations of motion are 

Rjlu = v a v (/ a )a - inv - iv M v,/i (153) 

R L = V^V v h uu - Dh - Ah 
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where is the covariant derivative associated with g. Note that h 
g^ v h^ v . We write explicitly 



R = (g^-h^)R^ (154) 



= R(o) + R L '■> 

where i?( ) is the Ricci scalar in the background. 

To linear order in the variation, the equations of motion then become 

[22] 

2A 2A 

5(Q^ + E^) = [1 + 2A(a + lp)\g%, + a[(D - —)G^ - —R L g^\ 

+(a + 2/3) [- V M V„ + g^D + Ag^]R L . (155) 

The most common gauge used in holography is radial axial gauge, h pfl = 
for n = (p,t,x,y). However, [22] used a covariant gauge 

V V = V„/i, (156) 

since in this gauge the equations of motion immediately simplify. 

Substituting (156) into the linearized tensors (153) gives 

9%, = \^^uh- ^□V + 7^V + £ A H«' ( 157 ) 
# L = -A/i 

which can then be substituted into (155), the variation of the equations of 
motion. Tracing over the result yields 

= 8{G txv + S^) = A[h - 2(a + 3P)Dh}. (158) 

Imposing the above constraint equation for the trace, we find that the vari- 
ation of the field equations is, 

o = 6{g^ + £^) = -|n 2 v-^(i + ^ + 8A/3^)nv 

A / 4Aa , A . 
+ 3 (1 + — + 8A/3J V (159) 

+ 3(a + 2 /3) (^-l M + A)v M V,/ i 
A (5a + 6(3 4A, „A A , 
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In [22] this equation was analysed only in the case of (a + 3/3) = 0, with a 
view to critical gravity, but we will not impose this constraint here. From 
(159) we wish to extract the equation of motion for h^, the traceless part 



of h^ u , where 



K v = h^ u) + ^g^v (160) 



This yields, provided that (3 / 0, 

This equation represents an inhomogeneous equation for the traceless part 
of the metric fluctuation. However one can rewrite the equation as a homo- 
geneous equation by defining a new traceless tensor ip^ u ) as 

Vv> = h <P») + AV <M V ^) /l ' ( 162 ) 

and choosing A such that the final term in (161) is zero. This value turns 
out to be 

A- 6{a + m (163) 
A -"3 + 8A(a + 3/3)' 

making the resulting equation of motion for -ip homogeneous 

0= ( D -i)( D -y- M2 )^>' (164) 

where 

M 2 = -2A-i-^. (165) 

a a 

Moreover it is easy to verify that ip^u) is transverse. 
In the case where a = 0, the analogue of (161) is 

(1 + 8/3A)(D - ^) V„) + (1 + 12/3A)V <M V„) = 0, (166) 

which can be rewritten as a homogeneous equation 
2A 

(□- T )^ M = 0; (167) 

6/3 „ „ , 
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Note that this equation is only second order. 

The interpretation of these equations is as follows. In the Einstein theory 
the only propagating mode is the traceless part of the metric, which couples 
to the dual stress energy tensor. In the theory with generic values of (a, f3) 
the trace of the metric is a propagating mode dual to a scalar operator Oh 
of dimension 

Ao » = f + /!+^T^ <168) 

whilst the equation of motion for the traceless part of the metric fluctuation 
is fourth order. One can write a basis for solutions of this fourth order 
equation as 

(□-|-M 2 )*a - »• 

with the propagating massless mode coupling to the dual stress tensor 

(2) 

and the new mode being associated with a spin two operator X of 

dimension 

Ax = ^ + ^v / 9 + M 2 . (170) 

In this section we will show explicitly how these modes are associated with 
the dual spin two operator and we will discuss how the variational problem 
is defined. 

Note that the special case in which the action on AdS is uncorrected, 
with the bulk term in the action reducing to Riemann squared, is obtained 
by choosing a = 47, j3 = —7, in which case 



^o h = o + a/7 + ^; ( m ) 




3 1 1 

A * " 5 + sf - ? 

We will discuss later when these operators are unitary, but let us already note 
here that for 7—7-0, which is indeed the case when the higher derivative 
corrections are small, either one or the other operator necessarily has a 
complex dimension and thus violates unitarity. 
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For special values of (a, /3) one has to look more carefully to obtain the 
spectrum. At a = the higher derivative term consists just of the Ricci 
scalar squared, and the only new propagating mode in the bulk is the trace 
of the metric fluctuation, dual to a scalar operator. At a + 3/3 = 0, when 
the bulk term reduces to Weyl squared, the metric trace is zero so there is 
no dual scalar operator but there is still a propagating spin two mode dual 
to a spin two operator. Whenever 

(1 + 2Aa + 8A/3) = 0, (172) 

the second spin two mode becomes massless, with the dual operator becom- 
ing the logarithmic partner of the stress energy tensor in the dual (L)CFT 
[47]. Note that this mode can become massless even when the trace is a 
propagating mode, with a + 3/3 = being an additional constraint used to 
remove the scalar operator. 

4.2 Derivation of equations of motion in general gauge 

In this subsection we derive an elegant form for the linearized equations 
of motion without imposing a gauge. Taking the trace of the equations of 
motion (41) to (43) one obtains 

(2a + 6P)DR - 12 - R = 0, (173) 

where we use the explicit value of the cosmological constant together with 
the Bianchi identity 

V%„ = \v v R. (174) 

Now letting r = (R+ 12) one obtains a diagonal equation of motion for r 

(2a + 6/3)Dr - r = 0. (175) 

Note that this equation did not rely on the linearized approximation and is 
exact. Defining 

Rfiv + 3g»v = + -rg^ (176) 

where is traceless, i.e. g^s^ = 0, the traceless part of the linearized 
equation of motion gives 

(«□ + (1 - 4a - 24/3)) = (a + 2/?)V (/1 V„ ) r, (177) 
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with the parentheses denoting the symmetric traceless combination. This 
equation can be diagonalized by defining 

= V - TV^V^r (178) 

with 

7 = 2 " + 3/? (179) 

7 3(1 -8a -24/3) 1 ' 



(<*□ + (1 - 4a - 240))V>L, = 0. (180) 



to give 

; , , n n _ i/, _ 9 i rw , 

The equations (175) and (180) represent second order equations for the 
linearized curvature tensor and hold in any gauge. In the covariant gauge 
V^/i^ = V u h used previously 

r = 3h; s^ v = —(□ + 2)h {H + iv^V^. (181) 

A complete basis for the solutions to (175) and (180) can be obtained by 
setting the metric fluctuation h^ u to be 

V = h lv + h*„, (182) 

with 

r{h T ) = 0; s^(h T ) = 0, (183) 

and (r{h x ),ij) r liU {h x )) are non-zero, satisfying (175) and (180). 

In understanding the holography dictionary it is useful to look at the 
asymptotic solutions for (175) and (180). Since r is simply a scalar field, of 
a specific mass, the general asymptotic solution to (175) is as usual 

r(p, x) = p 3 ~ A °r ( ro {x) + p 2 r 2 {x) + •••) + p A ° r (r 2 A-s(aO + ■ • • )> (184) 



where Ao r = § + yf + 2a+6/3 • H ere r o( x ) ac ts as the source for the dual 
operator, with r 2 A-3(a;) being the normalisable mode, and all other terms 
in the expansion being fixed by the field equation. 

Equation (180) is an equation for a massive spin two field of a given 
mass. Such fields are considered less frequently in holography (they were 
first analysed in [55]) but one can analyse the general asymptotic solutions 
to the field equations as follows. The independent solutions are 

iP {pp) (p,k) = p d - A (f(x) + ---)+p A (f(x) + ---) 

1>{i P )(p,k) = p d - A -\B l (x) + ---) + p A ~\B l (x) + ...) (185) 

%>(p,fc) = p d - A - 2 {X lJ {x) + ---) + p A - 2 (X lJ {x) + •••), 
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where 

A = i + ^d? + AM 2 (186) 

with d = 3 in this case and M 2 given in (165). The fields without tildes de- 
note the non-normalizable modes and those with tildes are the normalizable 
modes. Only the transverse traceless part of X^ and X^ are independent 
data, however, since the field equations imply 

X\ = Xf = 0- (187) 
Bi = -^Wjii 
Bi = - — _ - &Xij; 

f = ^—r&Bi\ f = -^d i B i . 

j 3- A A 

Thus the defining data for the spin two field indeed corresponds to a trans- 
verse traceless spin two operator in the dual field theory. 

The new defining boundary data is r^(x) in (184) and X^ in (185), 
namely the near boundary behaviour of the scalar curvature and the (trace 
adjusted) Ricci tensor. One can obtain a geometric interpretation of these 
boundary conditions as follows. As commented earlier, the most natural 
gauge for holography is the radial axial gauge in which the metric pertur- 
bations satisfy h pfl = and 

ds 2 = ^ + ^(Vij + H^dxUxi. (188) 

In this gauge the linearized Ricci scalar r[h] defined above is given by 

r[H] = p 2 R[H] - p 2 tr(H") + Spti(H'), (189) 

with a prime denoting a radial derivative and Rij being the linearized cur- 
vature of Hij , namely 

Rij = \ (d k d 3 H ik + d k d % R ik - didjXx{H) - 8 k d k H l3 ) ■ 
R = tfdjHij - □tr(.ff). (190) 

From the equation for the linearized Ricci scalar we note that the leading 
asymptotic behavior of the metric perturbation corresponding to the prop- 
agating scalar mode satisfies 

tlfl - (3-A)(l + A) ^ A °"-'°' + ^ <191) 
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We can also express this condition in a more geometric way, in terms of the 
extrinsic curvature of the hypersurface with induced metric r y fiu , as 



C K 



-4r 



(O)xP' 



3-A Gr 



(192) 



Therefore the new data r(o) (x) supplied corresponds to specifying the bound- 
ary condition for the trace of the normal derivative of the extrinsic curvature. 



4.3 Two point functions 

To extract the two point functions of the scalar and spin two operators, 
the field equations alone do not suffice: one needs to compute the onshell 
renormalized action. This is a non-trivial issue, as even when the bulk field 
has a mass such that the dual operator would be unitary, the corresponding 
two point function of that operator is not guaranteed to be positive. In 
other words, the sickness of the higher derivative theory can manifest itself 
in negative norms. 

A useful trick for obtaining the two point functions is the following, 
borrowed from the three dimensional discussions in [56] . Let us first rewrite 
the bulk terms in the action as 



1 

2^2 



d 4 x 



-9 



R - 2A + (f3 + j)R 2 + aS^S^ 



(193) 



where S^ u is the traceless part of the Ricci tensor. We next introduce a 
scalar auxiliary field <p and a traceless spin two auxiliary field ^ v and write 
the action as 



1 

2^2 



d 4 x 



-9 



a, 



2^1 



d 4 x 



R-2K+{l5+j){2R<j)- <p 2 ) 



(194) 



Eliminating the auxiliary fields using their equations of motion gives the 
previous action. For this action, the boundary term needed for a well-defined 
variational problem is 



Ighy = 2 / 



-7 



K 



a, 



1 + (2/3 + -)</> 



+ aK, 



(195) 



Note that the problems in setting up a variational problem have been solved 
here, by the introduction of the auxiliary fields. A similar approach to deal- 
ing with the variational problem in higher derivative theories was discussed 
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in [34]. The action with auxiliary fields admits Einstein manifolds as solu- 
tions, in which 

9nv = (f> = R = 4A; (p^ = S^ u = 0. (196) 

The action of course also admits other solutions, but in this section we are 
interested in the spectrum around an Einstein solution. For such solutions, 
the boundary counterterms given previously in (47) renormalise the onshell 
action. Note that, as previously anticipated, when one looks at the leading 
order Einstein solutions, the boundary conditions for the auxiliary fields do 
not involve non-trivial data (i.e. unlike the metric boundary condition, the 
boundary data for the auxiliary fields is not specified by arbitrary scalars or 
tensors) and indeed this remains true when evaluating corrections on such 
solutions. When we compute the spectrum below, however, we find that 
there is indeed non-trivial data required for the auxiliary fields, which is 
expressed in terms of arbitrary scalars and tensors. 

Let us now consider perturbations around such an Einstein solution g^ v 
of the equations of motion, i.e. we let 

Qixv = g^u + V/l ^ = 4A + 5(f); (p^ = 6(p^. (197) 

The boundary data for 5(f> and 5(f>^ v specify the defining data for dual scalar 
and tensor operators, respectively. 

Let us begin with the a = case. To quadratic order in the fluctuations 
one obtains the following for the bulk terms in the action 

51 = ~2^I d ' X ^ ^ ( G ^ [h] ) (198) 
+ 4 / d A xy^h^ (V^VM - U5<Pg^) 

K J 

+4 / d 4 x ^g 5(f) (R[h] -Ah- 5(f)) . 

where \i = (1 + 8/3 A), h = g^ u h^ v and the linearisation of the Einstein 
equation is 

QiM = R ^[h] - AV " 2 R ^ 9 ^ + \ Ah 9w ( 199 ) 
and the linearised Ricci tensor is given by 

R»v [h] = - (VV^V + vv^v " V„ V„/i - □ V) , (200) 



46 



with R[h] = g^R^u{h] being 

R[h] = V^V U V - Dh. (201) 
The action can be diagonalised with the field redefinition 

V = V + C(<¥)^ (202) 

and letting 

C = -2/3//U- (203) 
The bulk action at the quadratic level then becomes 

51 = dAx ^9h" U {^M) ( 204 ) 

-4 I d A x^^5<j)\-Q-U54> + 54> . 

k J La* 

where we have used 

R{6<f> 9lxl/ ) = -3DS<f>. (205) 

The equations of motion resulting from this action describe the graviton 
together with the scalar field, and agree with those found in the previous 
sections. 

Having obtained the action, it is now straightforward to extract the 
two point functions of the dual operators. To do this we need to keep 
careful track of the boundary terms, include those which arise in the field 
redefinitions. In working out these terms it is convenient to fix a gauge for 
the metric perturbation h^ u , the holographic radial axial gauge in which 

h pil = 0. (206) 

Thus the perturbed metric may be written as 



ds 2 = ^ + \ (rHj + Hij) d^dxi (207) 



where Hy = p 2 hij. Evaluating all boundary terms involving Hy, includ- 
ing those from the counterterms needed to renormalise the action for the 
background solution, one obtains 

/onshdi = -£2 / <$ x \ ( Hiid P H a + 1Hd p H) , (208) 
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where H = rf^Hij. These terms can be processed using the Fefferman- 
Graham expansion, namely 

Hij = H( )ij + p 3 H( 3 ) i:j H (209) 

with #(3)y being traceless and transverse. Thus 

/onshell = f (PxH'^H^. (210) 

This action is clearly finite, without the need for additional counterterms, 
as expected as the counterterms should take of all divergences of Einstein 
solutions of the field equations. Moreover, recalling that 

(T i7 -) = (an) 

we recover the formula 

( T ij) = J^2 H (3)iv ( 212 ) 

which is the linearisation of the renormalised stress tensor given earlier. 
Relative to Einstein gravity, this formula is shifted by a factor of fx which 
in turn implies that the two point function for the stress energy tensor 
will be shifted by factor of \i relative to the Einstein case [45]. In this 
theory the ratio rj/s is unchanged by the higher order correction. This was 
already apparent on general grounds, since the correction evaluated on an 
Einstein solution can be removed by field redefinitions. Here the derivation 
is somewhat non-trivial as both quantities are shifted by the factor of ji: the 
Wald entropy was computed earlier, and rj is obtained from the two point 
function of the stress energy tensor, which according to the formula above 
will only be shifted by n relative to the Einstein case. 

What remains is to collect together all of the terms involving the scalar 
field. These give 

/onshell = ~^2 J 3^)5^. (213) 

This is the action in Lorentzian signature. The corresponding action in 
Euclidean signature is 



Onshell = ^ / ' dEPd^SWto (214) 

liri J 



13/3 

[IK 

flK J p 
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Recalling that the asymptotic expansion of such a scalar field dual to an 
operator of dimension A is 

5<p = P d - A (5<f> id - A) + ...)+ P A (Ha + ■■■) (215) 

we see that this part of the onshell action still has divergences as p — > 0. 
This was indeed to be expected, as the counterterms computed earlier were 
for Einstein solutions of the field equations only. 

The holographic renormalization required for such a scalar field is already 
known: if the onshell (non-renormalized) Euclidean action for a free scalar 
field is 

iJL*& = -\J<Er<pVw (216) 

then the renormalised two point function of the operator of dimension A 
dual to the field </? is [57] 

= ffi r ^-$) * {£) s ■ (217) 

where we denote by 1Z the renormalised quantity. Comparing with our case 
we obtain the following for the norm of the two point function of the operator 
dual to 5(j): 

(0,(x)0,m = - 1 j^c A n^y (218) 

The norm is never positive and recall that the operator also has complex 
dimension for negative (3. It would be interesting to carry out a similar 
analysis for the spin two operator, setting a ^ 0, to find for which values of 
a the norm of the dual operator is non-positive. 

4.4 Spectra for other curvature corrections 

A similar analysis can be carried out for the spectrum around AdS induced 
by other curvature corrections. For any given curvature invariant one might 
anticipate that generically the operator associated with the higher derivative 
term is either non-unitary or has non-positive norm. There are certain 
exceptions to this generic case, however. 

If a curvature invariant which is built out of the Weyl tensor is added to 
the action, then the equations of motion linearised around AdS are necessar- 
ily unchanged since the Weyl tensor vanishes identically on the background. 
More precisely, any curvature at least cubic in the Weyl tensor implies that 
all contributions to the linearised field equations are at least linear in the 
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Weyl tensor of the background, which vanishes for AdS. Therefore the Weyl 
terms do not change the spectrum of operators in the dual CFT, although 
they can modify the correlation functions of these operators. This fits with 
the observation that the Weyl terms on asymptotically locally AdS space- 
times fall off sufficiently fast at infinity that the variational problem is un- 
changed from the Einstein case. Put differently, one needs no additional new 
boundary conditions and therefore there are no new associated propagating 
modes and corresponding dual operators. 

If one adds several different curvature invariants to the action, the diag- 
onalization of the linearised field equations becomes more complicated, as 
we saw in the case of curvature squared corrections. For each new boundary 
condition there is an associated new dual operator. The dimensions and 
norms of the dual operators are obtained non-trivially from diagonalising 
the field equations and manipulating the onshell action. 

From a top-down perspective, the leading higher derivative terms in the 
four-dimensional action must be consistent with unitarity. This implies that 
they must give rise to a linearised spectrum around AdS which is consistent 
with dual operators of real conformal dimension and positive norm. We have 
shown that individual terms such as R 2 are not consistent with unitarity, but 
we also noted that in reducing a higher dimensional curvature invariant any 
such curvature squared term always appears with fourth order curvature 
terms and a shift of the cosmological constant. Moreover, one needs to 
include all higher dimensional curvature invariants to respect unitarity at 
the required order, supersymmetry and so on. 

Finally let us consider how the dimensions of the dual operators relate to 
the parameters of the dual CFT. The dimensions of the operators for the case 
of Riemann squared were given in (171). From the discussion around (37) 
we note that when such a term arises from a reduction of eleven dimensions 
the coupling constant 7 would be the ratio of the term descending from 
R A to the leading order Einstein term. In other words, for the case of 
an S 7 reduction 7 would be of order 1/N and in ABJM it would be of 
order l/(kN) = 1/N'. The dimensions of the operators scale as l/y/^y, i.e. 
iV 1 / 2 and (iV') 1 / 2 = k A 1 / 2 respectively. Riemann squared on its own is not 
unitary, so one of the operators always has complex dimensions, but the 
combination with other terms arising from top down would give operators 
whose dimensions scale similarly. If no such operators exist in the dual CFT, 
then the net effect of the reduction of the leading top down terms must be 
trivial at the linearized level. 
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5 Conclusions 



In this paper we showed that the variational problem is not in general well- 
posed in higher derivative gravity theories without specifying additional data 
to the boundary metric. When the higher derivative terms are treated per- 
turbatively around the leading order Einstein solution, the higher derivative 
equations always become inhomogeneous second order equations, for which 
the variational problem is well-posed with only a boundary condition for 
the metric. However, in analyzing the spectrum around the corrected back- 
ground, the linearized equations of motion are generically higher order and 
do indeed require additional boundary conditions. In the context of hologra- 
phy these additional boundary conditions correspond to data for operators in 
the dual conformal field theory. For the curvature invariants we analyzed the 
operators are non-unitary since their conformal dimensions are generically 
complex and their norms are non-positive definite. From a top down per- 
spective, the reduction of any given higher dimensional curvature invariant 
results in a lower dimensional action involving several curvature invariants 
of different derivative order. When the lower dimensional curvature invari- 
ants are combined, the resulting spectrum must be unitary and thus either 
the dual operators must have real dimensions and positive norms or (per- 
haps more likely) the resulting lower dimensional linearized field equations 
remain second order with no new operators arising. 

Even when the new operators induced by the higher derivative terms 
are non-unitary, one might try to look for a unitary subsector of the theory, 
by switching off these operators. This is indeed the perspective of [21, 48], 
whose boundary conditions effectively set to zero the sources for the irrel- 
evant operators associated with the higher derivative terms. As explained 
in detail in [24, 25], however, switching off such sources does not switch 
off expectation values for such operators. Moreover, even if one restricts 
to a subsector of the theory in which the irrelevant operators are neither 
sourced nor acquire expectation values, the theory itself is non-unitary. In 
particular, the extra fields do contribute in computation of stress energy 
tensor correlation functions and there is no guarantee that the latter would 
be unitary, nor is it immediately apparent that the additional operators can 
always be decoupled from the stress energy tensor. 

Many interesting questions deserve further study. By comparison with 
dual field theory results, one could, at least in the highly super symmetric 
examples of AB JM and ABJ models, restrict what curvature invariants can 
arise in the effective four-dimensional action. It would also be interesting 
to work out the spectrum for the reduction of a top down curvature in- 
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variant, i.e. putting together curvature invariants of different order in four 
dimensions. In this work we have found that the simplest representative 
correction would involve the Weyl tensor at orders higher than two: such 
corrections are in some ways analogous to the Gauss-Bonnet examples in 
higher dimensions, in that the black holes are corrected but there are no 
new operators induced in the spectrum. One of the initial motivations for 
looking at higher derivative terms in AdS^ was to explore subleading effects 
in applied holography and the Weyl solution would be a natural candidate 
for such investigations. 

One approach to holographic cosmology exploits the domain wall cosmol- 
ogy correspondence [58] to obtain a field theoretic description of cosmologies 
in one higher bulk dimension [59]. Since the primary focus is naturally on 
four dimensional cosmologies, the results obtained here would be relevant in 
discussing higher derivative effects in holographic cosmology. In particular 
it would be interesting to understand whether corrections which give rise 
to non-unitary effects on the AdS side are automatically excluded from be- 
ing physical on the cosmological side, and whether the irrelevant operators 
associated with the higher derivative terms could actually be useful on the 
cosmology side in, for example, exiting from the inflationary era. 
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